THEORY OF THE PHASE TRANSITION 
FROM A DISORDERED CUBIC CRYSTAL 

TO A GLASS 

J. M. Yanezt, M. I. Molina^ and D. C. Mattis* 



^Facultad de Ciencias, Departamento de Ffsica, Universidad de Chile 

Casilla 653, Santiago, Chile. 
^Department of Physics, JFB Bldg. #201, 115 S. 1400 E 
Salt Lake City, UT 84112-0830, U.S.A. 

Abstract 

We calculate thermodynamic properties of a disordered model insulator, start- 
ing from the ideal simple-cubic lattice (g = 0) and increasing the disorder 
parameter g to ^> 1/2. As in the earlier Einstein- and Debye- approximations, 
the ground state energy is discontinuous at g c = 1/2. For g < g c the low-T 
heat-capacity C ~ T 3 whereas for g > g c , C ~ T. The van Hove singularities 
disappear at any finite magnitude g of the disorder. For g > 1/2 we discover 
novel fixed points in the self-energy and spectral density of this model glass. 
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The anomalous thermal properties observed in dielectric glasses p], 0, has 
prompted numerous theoretical and experimental studies on the subject [[|. 
Arguably, the best well-known model consists of the phenomenological two- 
level states (TLS) proposed by Anderson^] some time ago. Although the 
physical nature of the TLS remains elusive the theory does predict linear 
specific heat at low-T and enhanced ultrasonic attenuation, both universal 
properties of glassy materials 0. An alternative microscopic theory entirely 
based on phonons was subsequently advanced by one of the authors and his 
collaborators 0. It relies for its results on the interplay between disorder, 
parametrized by a dimensionless g, and anharmonicity, parametrized by a 
dimensionless £ (the Griineisen parameter) that is ultimately taken to 0. The 
solutions of two "toy" models, the one based on the unperturbed Einstein 
model of the phonon spectrum in the "crystalline" phase and the other, on 
the Debye model 0, [/J, both exhibited the following features: at a critical 
threshold g c the system undergoes a transition from a "disordered crystal" 
phase with partial long-range order (LRO) to a "glassy" phase characterized 
by a linear specific heat and a divergent Debye- Waller exponent-i.e. by a total 
loss of LRO0. 

Here, for the first time, we have extended this type of model disorder cum 
anharmonicity into the study of a bone-fide crystal and have obtained closed- 
form solutions in the simple-cubic lattice. This approach is the most realistic 
so far, given that initially the phonon density-of-states of the ideal crystal 
exhibits the van Hove singularities (vHS) thought to be characteristic of a 
crystal with LRO. The results are also much richer than in the toy models. In 
Fig. 1 we see that at low T our calculated specific heat is in better agreement 
with experiment than is the specific heat from molecular dynamics. 

Principally our findings are as follows: In the disordered crystalline phase, 
where the vHS might be thought to weaken with increasing g and ultimately 
disappear, instead they soften and disappear immediately at any finite g > 0, 
despite the persistence of LRO and of a finite Debye- Waller exponent up until 
g > gc- There is an exothermic first-order phase transition as g is increased 
above g c = 0.5. Nevertheless the elastic properties in the glassy phase are 
intimately related to those of the crystal, as we shall show. Additionally, in 
the glassy phase g > g c , we discovered that the spectral density function p 
exhibits fixed points: frequencies at which the spectral density is independent 
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of g. This feature bring us close to the ultimate goal of being able to express a 
microscopic law of corresponding states for glassy materials. The specific heat 
exhibits a Debye T 3 law throughout the crystalline phase, up to g c . The exact 
coefficient B(g) is given in closed form below. Beyond g c the low-T specific 
heat is linear in T, with coefficient A(T) also given in closed form below. 

THE HAMILTONIAN. In an effort to make this paper self-contained 
we review here the main features of the model. The starting point is the 
Hamiltonian of the reference crystal 

#o = $>* (4a* + 5) (fr = 1 )> (!) 

where k contains the wavevector and polarization of the phonons. The opera- 
tor a k (a&) are the familiar creation (destruction) phonon operators. Disorder 
is introduced in the model via the rather general phonon scattering term 

H x = j £ M(k, k') ^JU^Q k Qy , (2) 
4 k,k' 

where Qk = a_ k + ak and the M(k, k') are random scattering amplitudes 
which depend on the location and type of disorder. We model M(k, k') by a 
set of random phases (g/yN) exp[i 6(k, k')], where the 9(k,k') = —9(k',k) 
are random and distributed uniformly in [0, 2n] and g measures the strength of 
the disorder (g 2 is proportional to the concentration of defects in the crystal). 
This choice is what allows the calculation of the exact spectrum of eigenvalues 
of the random matrix||. Anharmonic effects are introduced next in order to 
stabilize the spectrum, since H + Hi possess an unstability threshold where 
the phonons become overcoupled^], [Tj. A simple quartic term suffices for this 
purpose, 

£-77 / ^k^k'QkQ-kQk'Q-k' , (3) 

where £ is the Gruneisen parameter, typically a small quantity, expressed 
here in dimensionless form. iV is the number of normal modes and Uq is a 
characteristic frequency of the reference crystal. Then a symmetry-breaking 
transformation — > Q\ + (/fc/^ 1 ^ 2 ) is performed in order to maintain all 
the u; 2 's positive. This shift, required for g > 1/2 to restore stability in the 
dynamics, produces new terms linear in Q's that have coefficients 0(1/ 1; 1 / 2 ) 
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but that arise from two distinct sources: from the quadratic terms (linear in 
the /'s) and from the quartic (Griineisen) terms (third-order in the /'s). These 
must be canceled if the individual atoms are all to be in positions of stable 
equilibrium. The relevant equations admit only a trivial solution f k = for 
g < 1/2, but for g > 1/2 the cancellation is effected by choosing f k to be an 
eigenfunction of the random scattering matrix. (To minimize the free energy 
one picks the eigenvalue belonging to the largest eigenvalue.) The ground state 
energy is then a discontinuous function of g at g c — 1/2. At this value of the 
disorder the original locations of the atoms are no longer positions of stable 
equilibrium. Each atom moves to new, random, position that is essentially 
uncorrelated with the motion of the other atoms. This picture is confirmed 
by a calculation of the Debye- Waller factor; we find it to be finite below g c 
but divergent (oc l/£) once g exceeds g c , thereby confirming that all long- 
range correlations are lost. The dynamical Hamiltonian|| |7j governing lattice 
vibrations is, however continuous, and independent of £ in the limit. It is: 



H = (a\a k + + V^k^k'Mik, k')Q k Q k > 

k v z/ 4 k,k' 

+\ (g - i) <? (g - ~) E <*kQkQ-k + O(0 + 

9-l){o^)-o((g-l) 2 /ty (4) 

where i? is the Heaviside step function. The origin of the individual powers in 
£ in H is as follows. 

The ^-independent correction to the elastic forces is the result of all contri- 
butions second-order in the /'s in the quartic (Griineisen) term. The singular, 
discontinuous, shift (oc l/£) in the ground state energy comes from the quartic 
term, when all four factors are /'s and from quadratic terms in which both 
factors are /'s. But, however large it may be, this shift in ground state energy 
does not affect the dynamics. 

Two remaining contributions to the transformation vanish in the lim £ — ► 0: 
terms with a single factor /'s in the quartic term yield contributions are 0(£ 1//2 ) 
and those that are quartic in Q's (with no /'s) yield the original quartic term 
°c£. 

The free energy of the model can be expressed as an integral over the 
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Po(u) = — w Im 

7T 



(6) 



coupling constant, 

F = k B T dw[po(u)+ Pl (u)+ P2 (u) ] log(2 sinh(/?cu/2)), /3 = l/k B T 
Jo 

(5) 

where poiyj) is the density of states of the reference crystal, 

r d 3 k 1 
Ifbz (27r) 3 uj 2 — Uj\ 

and the densities of modes p\{oS) and p%{tjS) are the contributions from disorder 
and anharmonicity[[7|: 

2 r 1 d\ f) 

Pi(f») = —z / ^^[-/(^A)^^)] (7) 

7T g z JO A 13 OUJ 

, , ^(o- 1/2) ri dA . , 9 r r , x . . 

and R(u,X), J (a;, A) are the real and imaginary parts of the the self energy 
Z(u,X) which obeys the transcendental equation || 

A) = ? - W 2 ( i + J + ^ A)) (9) 
with = (2A# - 1) $(A# - 1/2). 

RESULTS FOR THE SIMPLE-CUBIC LATTICE. The integral 
(||) can be performed analytically if the dispersion relation are io\ = Uq(3 — 
cos(fci) — cos(k 2 ) — cos(/c 3 )). After inserting this dispersion into Eq.(|9]) and 
taking the continuum limit, the equation for the self energy Z(u, A) can be 
written implicitly as, 



W [3, r(w)] = (l - 3^|(1 + cj> + Z(u))\ (l 
where the argument t(uj) is given by, 
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■(1 + + zh) , (io) 



r(lA 3(l + + ZM) 

1 ' 3(l + a^))-W 1 } 

and W[3, z] is the generalized Watson integral evaluated by Joyce ||, who ex- 
pressed it compactly in terms of complete elliptic integrals of the first kind|T0|: 



Tr[;J _ :] j'2\ 2 \J^^ K{k+)K{k _ ) (12) 

\7T/ 1 — X\ 
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where x x = (1/2) + (z 2 /6) - (1/2) V / T^^^1 - (z 2 /9),x 2 = x x /{x x - 1), 
k% = (1/2) ± (l/4)x 2y /4 -x 2 - (1/4) (2 - x 2 )y/l-x 2 . Figures 1 and 2 show 
the real and imaginary parts of Z(u,g). The imaginary part I(u,g) shows 
the not unexpected systematic broadening of the bandwidth with increasing 
disorder. R(u,g) exhibits a change in curvature from convex to concave at 
u = as g is increased from below to above g c = 0.5. Correspondingly, the 
slope of I(w,g) vanishes at u = for all g < 1/2 but acquires a finite value 
for g > 1/2. 

Z(u,g) can be expanded at both the low and high frequencies. It is also 
seen to possess a low-frequency fixed point in the glassy phase. From Eq.(|), 
for A > 1 and g > 1/2 we deduce a sort of duality relation, 



connecting high frequencies, large disorder, to lower frequencies and smaller 
disorder. Z is inserted into the global density-of-states p = po + pi + p 2 and the 
indicated integrals performed, with the results shown in Fig. 3. For g > 1/2, 
two fixed points in p are discerned: one just below the characteristic frequency 
ujq and the other just below 3ujq. A new high-frequency "tail" grows beyond 
the second fixed point, indicating an accumulation of spectral density at the 
highest frequancy as well as at the lowest. 

It is also clear that the two van Hove cusps at to/cuo = \/2 and 2 are sharp 
only at g — 0, and that p losses its van Hove cusps and becomes analytic in 
uj for any nonzero value of the disorder, i.e. already in the disordered crystal 
and not just in the glassy phase! The theory works extraordinarily well for g 
from zero into the glassy phase at 0.5 + e, with almost every noncrystalline 
material being excellently characterized at low T by some e < 0.1, and it 
remains trustworthy in every detail until g exceeds 1, when the mean- field 
approximation to the quartic terms fails at high frequencies. 

THE SPECIFIC HEAT. With a knowledge of p the model's thermody- 
namic properties, including specific heat, entropy and the other such functions 
are readily obtained. For g < 1/2, we find at low T, 



Z(u,\g) = A Z(uj/VX,g) 



(13) 



c = 




(14) 
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where B(g) = (4n 2 V2/5) ( 1 + 7 ( [ (1/2) + ^(1/4) - ^ 2 ]- 3 / 2 - 1 ) ) and 
7 ss 1.43. For g > 1/2, 



where A(<?) = tt(2/3) 5 / 2 v / W[3, 1](1 - (2<?)- 3 / 2 ) « 1.404(1 - (2<?)- 3 / 2 ). 



The high-T specific heat is reasonably universal, always tending smoothly 
to the Dulong- Petit limit because / °° duo p(uo) = 3. Figure 4 illustrates the sep- 
arate contributions Co, C\ and C2 from the crystal, disorder and anharmonicity 
respectively, at g — 0.7, over the interval < T/T < 0.6. 

SPEED OF SOUND. From our calculations using the present model 
and on the earlier Debye model, we find that, at constant density, the long 
wavelength speed of sound s = (du/dk) k ^ Q , when expressed as a ratio to the 
speed of sound in the perfect crystal, v(g) = s(g)/s(0), is a universal function 
of g having a cusp minimum at g — 1/2. Explicitly, we obtain: 



identical to Fig. 2 in Molina and Mattis [6] for the Debye model. The predicted 
drop of some 28% as g increased from to 1/2 + e should be confronted 
with experiment, after corrections for any changes in density with increasing 
disorder are taken into account. 

We intend to publish the mathematical details of the various calculations 
elsewhere, but in the meantime the interested reader may request them from 
the corresponding author, M.I.M. 
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for g < 1 /2 
for g > 1/2 



7 



References 



[1] R. C. Zeller and R. O. Pohl, Phys. Rev. B 4, 2029 (1971). 

[2] R. O. Pohl, in W. A. Phillips (Ed.), Topics in Current Physics, Vol. 24, p. 
26 ( Springer- Verlag, 1981). 

[3] R. Zallen, The Physics of Amorphous Solids (Wiley, New York, 1983); 
W. A. Phillips, (Ed.), Amorphous Solids. Low Temperature Properties 
(Springer, Berlin 1981). 

[4] P.W. Anderson, B. Halperin and C. Varma, Philos. Mag. 25, 1468 (1972); 
W. A. Phillips, J. Low Temp. Phys. 7, 351 (1972). 

[5] S. Hunklinger, W. Arnold, S. Stein, R. Nova and D. Dransfeld, Phys. Lett. 
A 42, 253 (1972). 

[6] Daniel C. Mattis, Phys. Lett. A 117, 297 (1986); Daniel C. Mattis, Phys. 
Lett. A 120, 349 (1987); M. I. Molina and D. C. Mattis, Phys. Lett. A 
159, 337 (1991); Mario Molina, Phys. Lett. A 161, 145 (1991). 

[7] Mario Molina, Ph.D. thesis, University of Utah, 1991. 

[8] See D. Mattis and R. Raghavan, Phys. Lett. 75 A, 313 (1980) and refer- 
ences therein. For present purposes we require only the largest eigenvalue, 
obtainable by simple random-walk theory. The actual density of random 
eigenvalues is irrelevant, for as soon as the largest eigenvalue exceeds a crit- 
ical value the normal mode spectrum becomes unstable and atoms move 
to their new positions of equilibrium. 

[9] G. S. Joyce, J. Phys. A: Gen. Phys. 5, L65 (1972); G. S. Joyce, Philos. 
Trans. A 273, 583 (1972). 

[10] Milton Abramowitz and Irene Stegun (Eds.), Handbook of mathematical 
functions, Dover Publications, New York, 1972. 



8 



Captions List 



Figure 1: Detailed comparison of specific heat of a-SiC>2 with the present 
theory and with computer simulation molecular dynamics (S. N. Taraskin and 
S. R. Elliott, Phys. Rev. B 56, 8605 (1997)). We observe excellent agreement 
of the present theory with experiment at the low -and high- temperature end, 
although our result does lack the "bump" at 10 K. 

Figure 2: Real and imaginary parts of the self-energy Z(u) as a function of 
frequency for several disorder concentrations, below and above critical (g = 
1/2). 

Figure 3: Total density of states p(u) = pa{oj) + pi(u) + P2^) as a func- 
tion of frequency for a wide range of disorder concentrations, proportional 
to g 2 . Arranged according to the height of the central maximum, we have 
= 0, 0.1, 0.2, ... down to g = 1. 

Figure 4: The three distinct contributions to the low-temperature specific 
heat from Eqs. (^|), (|7|) and (§) in the glassy phase (g = 0.7). 
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